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ABSTRACT 



Breakup temperatures were deduced from double ratios of isotope yields for target 
spectators produced in the reaction 197 Au + 197 Au at 1000 MeV per nucleon. Pairs of 
3 ' 4 He and 6 ' 7 Li isotopes and pairs of 3,4 He and H isotopes (p, d and d, t) yield consistent 
temperatures after feeding corrections, based on the quantum statistical model, are applied. 
The temperatures rise with decreasing impact parameter from 4 MeV for peripheral to 
about 10 MeV for the most central collisions. 

The good agreement with the breakup temperatures measured previously for projec- 
tile spectators at an incident energy of 600 MeV per nucleon confirms the universality 
established for the spectator decay at relativistic bombarding energies. The measured 
temperatures also agree with the breakup temperatures predicted by the statistical multi- 
fragmentation model. For these calculations a relation between the initial excitation energy 
and mass was derived which gives good simultaneous agreement for the fragment charge 
correlations. 

The energy spectra of light charged particles, measured at 9^ = 150°, exhibit Maxwellian 
shapes with inverse slope parameters much higher than the breakup temperatures. The 
statistical multifragmentation model, because Coulomb repulsion and sequential decay pro- 
cesses are included, yields light-particle spectra with inverse slope parameters higher than 
the breakup temperatures but considerably below the measured values. The systematic 
behavior of the differences suggests that they are caused by light-charged-particle emission 
prior to the final breakup stage. 



Keywords: 197 Au projectiles and targets, E/A = 600 and 1000 MeV; measured fragment 
cross sections, isotopic yield ratios; deduced breakup temperatures, pre-breakup emission; 
analysis using quantum statistical and statistical multifragmentation models. 

PACS numbers: 25.70.Mn, 25.70.Pq, 25.75.-q 
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1 Introduction 



Heavy ion reactions at relativistic energies offer a wide range of possibilities to study the 
multi-fragment decay of highly excited nuclei [1-8] . In collisions of heavy nuclei at incident 
energies exceeding values of about 100 MeV per nucleon , highly excited and equilibrated 
spectator systems are formed which decay by multifragmentation M in good agreement 



with statistical predictions [10, 11]. The analysis of the kinetic energies of the decay 



products has not revealed significant flow effects 0, PS], O. Therefore, the spectator nuclei 



which are produced over wide ranges of excitation energy and mass in these reactions, 
are well suited for the investigation of highly excited nuclear systems in thermodynamical 
equilibrium. 

From the simultaneous measurement of the temperature and the excitation energy for 
excited projectile spectators in 197 Au + 197 Au collisions at 600 MeV per nucleon, a caloric 



curve of nuclei has recently been obtained | 13fl . For the temperature determination the 



method suggested by Albergo et al. has been used which is based on the assumption of 
chemical equilibrium and requires the measurement of double ratios of isotopic yields fl4 



The obtained temperatures, plotted against the measured excitation energy, resulted in a 
caloric curve with the characteristic behavior reminiscent of first-order phase transitions in 
macroscopic systems. The 'liquid' and the 'vapor' regimes where the temperature rises with 
increasing excitation energy are separated by a region of nearly constant temperature T rs 
5 MeV over which the multiplicity of the fragmentation products increases. These results 
and, in particular, the apparent rise of the breakup temperature at excitation energies 
exceeding 10 MeV per nucleon have initiated a widespread discussion which addresses 
both methodical aspects and questions of interpretation [15-22] . 

The qualitative shape of the caloric curve of nuclei has been predicted long ago on the 



basis of the statistical multifragmentation model [2^] . The same model, more recently, has 
been shown to rather accurately describe the charge correlations measured for the reaction 
197 Au on Cu at E/A = 600 MeV, including their dispersions around the mean behavior fPT |. 
Naturally, the question arises whether a comparably quantitative level of accuracy can be 
reached for the reproduction of the measured caloric curve. A statistical description of the 
spectator decay will only be consistent if the model parameters, including the temperature, 
are in the range given by the experiment. The breakup temperature should also exhibit the 
invariance with respect to the entrance channel that has been found for the fragmentation 
patterns. Their universal dependence on Zb oun d (Zbound scaling) is a prominent and well 
established feature of the spectator decay where Zbound, defined as the sum of the atomic 
numbers Z$ of all projectile fragments with Zi > 2, is a quantity closely correlated with 
the excitation energy imparted to the primary spectator system |§. 

In this work, we present results of a new measurement of the breakup temperature in the 
197 Au + 197 Au reaction, and we address these open questions associated with the statistical 
interpretation of multi-fragment decays of excited spectators. High-resolution telescopes 
were used to measure isotopically resolved yields of light charged particles and intermediate- 
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mass fragments emitted by target spectators produced at a bombarding energy of 1000 MeV 
per nucleon. Methodically we relied upon the observed isotropy of the spectator decay || 
and replaced the measurement of solid-angle-integrated yields by that of differential yields 
at selected angles. Temperatures were deduced from double yield ratios of H, He, and Li 



isotopes whereby feeding corrections, obtained from the quantum statistical model [[24| [25* 
were applied. 

The deduced breakup temperatures are compared to those measured previously with the 
ALADIN spectrometer for the 197 Au + 197 Au system at 600 MeV per nucleon |26| and to 
predictions of a recent version of the statistical multifragmentation model |27| . The input 
parameters for these calculations were derived from the requirement of a good simultaneous 
reproduction of the fragment charge correlations. While these comparisons turn out to be 
rather favourable, it is evident from the measured energy spectra that emission prior to 
the multi-fragment breakup contributes significantly to the observed light-particle yields. 



2 Experimental method 

The experiment was performed at the ALADIN spectrometer of the GSI facility ]9 
Beams of 197 Au of 1000 MeV per nucleon incident energy were provided by the heavy- 
ion synchrotron SIS and directed onto targets of 25-mg/cm 2 areal thickness. The present 
data were taken as part of a larger experiment which incorporated three multi-detector 
hodoscopes for correlated particle detection, built at INFN-Catania, GSI, and Michigan 
State University and consisting of a total of 216 Si-CsI(Tl) telescopes [^9|. |3(J. 

A set of seven telescopes, each consisting of three Si detectors with thickness 50, 300, 
and 1000 /im and of a 4-cm long CsI(Tl) scintillator with photodiode readout, were used 
to measure the isotopic yields of light charged particles and fragments. Four telescopes 
were placed in the forward hemisphere while three telescopes were placed at d\ a b = 110°, 
130°, and 150° for detecting the products of the target-spectator decay. Each telescope 
subtended a solid angle of 7.0 msr. Permanent magnets were placed next to the entrance 
collimator of each telescope in order to deflect 5 electrons emerging from the target. 

For a global characterization of the reaction and impact parameter selection, the quan- 
tity Z bound of the coincident projectile decay was measured with the time-of-flight (TOF) 
wall of the ALADIN spectrometer. Because of the symmetry of the collision system, the 
mean values of Zb OU nd, for a given event class, should be the same for the target and the 
projectile spectators. Only very small differences may arise from the finite dispersion of 
the relation between Zbound and the impact parameter |5I|. 



The energy calibration of the silicon detectors was obtained from the calculated punch- 
through energies for several isotopes. Radioactive sources of a particles were also used. 
For the CsI(Tl) detectors, energy calibration was achieved by using the calculated punch- 
through energies of protons, deuterons, and tritons, the calculated energies of particles 
whose energy loss was measured with the preceding Si detector, and theoretical light- 
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output curves |32|, |33|| . 



The quality of the obtained particle identification is illustrated in fig. 1 for the case of 
the telescope placed at 9i a b = 150°. Isotopes in the range from hydrogen to carbon were 
satisfactorily resolved. Isotopic yields were determined by Gaussian fitting with background 
subtraction. The yields of particles not stopped in the CsI(Tl) detectors was very low at 
these backward angles and did not affect the deduced isotope ratios. The contamination 
of lithium yields by two- 4 He events (decay of 8 Be) was estimated to be 2% and ignored. 
However, because of the presumably different shapes of 3 He and 4 He spectra at low energies 
(cf. refs. E2, [351]), a correction was required in order to compensate for the effect of 



the detection threshold E 12 MeV for helium, resulting from triggering with the 300-/im 
detector. The extrapolation to zero threshold energy was based on Maxwellian spectra 
shapes fitted to the measured parts of the spectra (fig. 2, full lines). The systematic 
uncertainty associated with this procedure was assessed by comparing to other methods 
of extrapolation. An upper limit of the 3 He/ 4 He yield ratio was obtained from a data set 
composed of events triggered by other detectors and extending down to the identification 
threshold for helium of E m 8 MeV. A lower limit was obtained by choosing different 
energy thresholds E > 12 MeV in the off-line analysis and by linearly extrapolating the 
resulting isotope ratios to zero threshold height. The linear extrapolation overestimates 
the threshold effect because the particle intensities should decrease at very small energies. 
For these extrapolations of the helium yield ratio, the fitted and background corrected 
yields as obtained from the identification spectra were used. 

For the hydrogen and lithium isotopes, extrapolations to zero threshold were not consid- 
ered necessary and the measured yield ratios above threshold were used. The experimental 
thresholds are sufficiently low for the hydrogens and, in the lithium case, a significant mass 
dependence of the energy spectra was neither expected nor found. This is illustrated in fig. 
3. The curve derived from fitting the spectrum of 7 Li ions describes very well also that of 
6 Li (left panel). The 6 Li/ 7 Li yield ratio does not change significantly as a function of the 
applied threshold energy (right panel). 



3 Experimental results and discussion 

3.1 Experimental temperatures 

Emission temperatures T were deduced from the double ratios R of yields Y{ mea- 
sured for pairs of neighboring H, He, or Li isotopes. Under the assumption of chemical 
equilibrium, they may be expressed as 

R = ^ = a ■ exp(((5 1 - B 2 ) - (B 3 - B A ))/T). (1) 

Y 3/ Y 4 

Here Bi denotes the binding energy of particle species i and the constant a contains the 
ground-state spins and mass numbers. For the ratios to be sufficiently sensitive to temper- 
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ature the double difference of the binding energies 



- B 2 ) - (B 3 - B A ) (2) 



should be larger than the typical temperature to be measured |[20|| . In this work we chose 
3 He/ 4 He as one of the two ratios (the difference in binding energy is 20.6 MeV) which 
was combined with either the lithium yield ratio 6 Li/ 7 Li or with the hydrogen yield ratios 
p/d or d/t. The set of 3 He, 4 He, 6 Li, and 7 Li isotopes is the one used previously for the 
determination of breakup temperatures of projectile spectators in 197 Au + 197 Au at 600 
MeV per nucleon ||13|| . Combinations involving p, d, or t, together with 3 He and 4 He, have 
the advantage of larger production cross sections, particularly in the 'vapor' regime where 
the heavier fragments are becoming rare. 

The four yield ratios used in this work are shown in fig. 4 as a function of Zbound- For 
3 He/ 4 He the systematic errors corresponding to the two limits of extrapolation, described 
in the previous section, are indicated by the brackets. The particular sensitivity of this 
ratio is evident; with decreasing Z bound , i.e. increasing centrality, it increases by about one 
order of magnitude while the other three ratios exhibit a weak and qualitatively similar 
dependence on Z bound . 

Solving eq. (1) with respect to T yields the following three expressions: 

T HeLi ,o = 13.3/ln(2.2 ^ Li (!; Li ) (3) 

* 3 He/ 1 4 He 

and 

T Hep d,o = 18.4/ ln(5.6 ^ H ^ H ) (4) 

"He/ "He 

and 

T H edt,o = 14.3/ln(1.6 ^ H ^ 3H ) (5) 

"He/ "He 

where the temperatures are given in units of MeV. The subscript is meant to indicate 
that these apparent temperatures, derived from the measured ground-state populations, 
may be affected by feeding of these populations from sequentially decaying excited states. 
The consequences of side feeding are presently a subject of considerable discussion and 
are quantitatively investigated by several groups with different models [36-39]. Here, the 
required corrections were calculated with the quantum statistical model which starts from 
chemical equilibrium at a given temperature, density, and neutron-to-proton (N/Z) ratio 
and which includes sequential decay p|, p5| . In fig. 5, the three apparent temperatures 



defined in eqs. (3-5) are shown as a function of the equilibrium temperature T input for 
the parameters N/Z = 1.49 (value of 197 Au) and density p = 0.3 • po (where po is the 
saturation density of nuclei). The relations between TH e Li,o or TH e dt,o and T input are almost 
linear and the corrections required in these two cases are practically identical, except at the 
highest temperatures. The linear approximation, indicated by the dotted line, corresponds 
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to the constant correction factor T = 1.2 ■ T adopted for T He Li in ref. (T^j. The figure 



also demonstrates that TH ep d,o is, apparently, more strongly affected by feeding effects at 
the higher temperatures. In this case, the double ratio includes the yield of protons which 
are likely to be produced in the decay of excited light fragments. Within the range of 
densities 0.1 < p/p < 0.5, the corrections required according to the quantum statistical 



model vary within about ±15% [26]. They virtually do not change with the N/Z ratio of 
the primary source. In the analysis, the corrections calculated for p/po — 0.3, as derived 
from the results shown in fig. 5, were applied. 

The side-feeding predictions obtained by other groups [36-39] suggest that the correc- 
tions are model dependent. In the region of low excitation energies, the correction factors 
required for T^Li range from 1.0 to 1.3, within the uncertainties of the model assumptions. 
However, at high excitation energies where the consequences of side feeding from higher 
lying states may become more important, the results differ considerably and depend on 
the amount of unbound states in the continuum that are considered 137|, |39| and on the 



assumed breakup density |26|, [38], |39|]. In the present work, we use the corrections obtained 



from the quantum statistical model in order to be consistent with previous analyses. We 
will also refer to the statistical multifragmentation model which predicts a very similar 
relation between T^Li and the equilibrium temperature (section 3.2). 

In fig. 6 the obtained temperatures Tn e Li, ^Hepd, and Xkedt are shown as a function 
of Zbound- They are based on the yield ratios measured with the telescope at the most 
backward angle Q\ a b = 150°. Simulations indicate that, at this angle, contributions from 
the midrapidity source should be small. Zbound was determined from the fragments of the 
projectile decay measured for the same event. The temperatures increase continuously 
with decreasing Zbound from T = 4 MeV for peripheral collisions to about 10 MeV for 
the most central collisions associated with the smallest Zb OU nd values. The range Zbound < 
20 corresponds to the high excitation energies at which the upbend of the temperature 
appears in the caloric curve [|T3|]. We note here that the cross section for Z bound < 2 is 
of the order of 130 mb in this reaction [0] so that the minimum impact parameter for 
Zbound > 2 (partitions with at least one fragment of Z > 2) is more than 2 fm in a sharp 
cutoff approximation. 

The results obtained with the three different double ratios agree rather well. Only the 
strong rise of T He Li at small Z bound is not equally followed by the other two temperatures. 
TneLi and Tfiepd track each other rather closely which is remarkable in view of the different 
feeding corrections (fig. 5). Only Tn e dt is systematically somewhat lower than the other 
two temperatures. Tn c dt has also been reported to be lower than T^Li in ref. [pTH , in 



agreement with the statistical-model calculations of ref. P5|, but has been found to be 



slightly higher than T He Li m ref- @ • It seems difficult to assess the precise nature of these 
deviations at the present time. 

A comparison of the TneLi temperatures with those derived in previous work for pro- 
jectile spectators in 197 Au + 197 Au collisions at 600 MeV per nucleon is given in fig. 7. In 
the case of the projectile decay, the isotopes were identified by tracking of their trajecto- 
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ries with the upgraded TP-MUSIC detector and subsequent momentum and time-of-flight 



analysis |TB|, The displayed data symbols represent the mean values of the range of 
systematic uncertainties associated with the two different experiments while the errors in- 
clude both statistical and systematic contributions. The projectile temperatures are the 
result of a new analysis of the original data and are somewhat higher, between 10% and 
20%, than those reported previously Their larger errors follow from a reassessment of the 
potential 4 He contamination of the 6 Li yield caused by Z misidentification. 

Within the errors, good agreement is observed for the results at 600 and 1000 MeV 
per nucleon. The expected invariance of the breakup temperature with the bombarding 
energy is thus confirmed. It is consistent with the Zbound scaling of the mean fragment 
multiplicities and charge correlations and supports the statistical interpretation of the 
multi-fragment decay of highly excited spectator nuclei 0. The breakup temperatures 
deduced by the EOS collaboration for 197 Au + C at 1 GeV per nucleon are also consistent 
with this conclusion [7J]. Within the range Zb oun d > 40 which is mainly populated in this 
reaction ||, they are in agreement with the present results for 197 Au + 197 Au, both in 
absolute magnitude and in their dependence on the impact parameter, and thus confirm 
the expected invariance with respect to the mass of the collision partner. 

3.2 Model calculations 



The calculations within the statistical multifragmentation model |27| were performed 
in order to test its consistency with respect to the statistical parameters and predicted 
charge partitions. Here one assumes that all observed particles come from the decay of 
one equilibrated source. The correlation of excitation energy and mass of the ensemble 
of excited spectator nuclei, required as input for the calculations, was chosen in the form 
shown in fig. 8. The neutron-to-proton ratio was taken to be that of 197 Au (N/Z = 1.49). 
The distribution of excitation energies at fixed spectator mass Aq had a Gaussian width, 
chosen in proportion to the square root of the mean excitation energy, and the distribution 
of masses A Q was adjusted in order to reproduce the measured cross section da / dZ hound . 

As a criterion for setting the excitation energy per nucleon, E x /A, we chose the ca- 
pability of the model to simultaneously describe the correlations of the mean multiplic- 
ity (Mimf) of intermediate-mass fragments (IMF's) and of the mean charge asymme- 
try (ai 2 ) with Zbound- The asymmetry a i2 of the two largest fragments is defined as 
a 12 = [Z max — Z 2 )/(Z max + Z 2 ), with the mean value to be calculated from all events 
with Z max > Z 2 > 2. The comparison, shown in fig. 9, was based on the data reported 
in ref. [§. In the region Z boU nd > 30, the mean excitation energy of the ensemble of 
spectator nuclei was found to be well constrained by the mean fragment multiplicity alone. 
At Zbound ~ 30 and below, the charge asymmetry was a necessary second constraint (cf. 



ref. ||40| ) while, at the lowest values of Z boU nd, neither the multiplicity nor the asymmetry 
provided rigid constraints on the excitation energy. These sensitivities are illustrated in 
fig. 9 where the dashed and dotted lines show the model results for E x /A chosen 15% 
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above and below the adopted values. 

The excitation energies that have resulted from this procedure are somewhat larger 
than those found previously in analyses [11, 41-43] of the earlier 197 Au on Cu data at 600 
MeV per nucleon fl~3| . The difference reflects the sensitivity to the fragment multiplicity 



and is caused by the slightly larger mean multiplicities that were obtained from the more 
recent experiments with improved acceptance |J. The excitation energies are still smaller 
than the experimental values obtained with the calorimetric method of summing up the 
kinetic energies of the product nuclei and their mass excess with respect to the ground 
state of the original spectator system M O . 



The calculations proceed such that, for a given mass, charge, and excitation energy, 
first the partition function is calculated and then the temperature T = (Tf) is obtained 
as the average over the ensemble of partitions. In the microcanonical approximation, the 
temperature Tf for a particular partition / is found from the energy balance 

E f (T f ,V)=E (6) 

where Ef is the energy of the system at Tf within the volume V, and Eq is the total 
available energy [M . At high excitation energies, when the system disassembles into many 
fragments, this temperature is very close to the grand canonical temperature. 

The solid line in fig. 7 represents the thermodynamical temperature T obtained in 
this way from the calculations. With decreasing Z bmin( i, h increases monotonically from 
about 5 to 9 MeV. Over a wide range of Zbound it remains close to T = 6 MeV which 
reflects the plateau predicted by the statistical multifragmentation model for the range of 



excitation energies 3 MeV < E x /A < 10 MeV p7 |. In model calculations performed for 
a fixed spectator mass, the plateau is associated with a strong and monotonic rise of the 
fragment multiplicities. Experimentally, due to the decrease of the spectator mass with 
increasing excitation energy, the production of intermediate-mass fragments passes through 
a maximum in the corresponding range of Zb oun d of about 20 to 60 (cf. figs. 8, 9). 

The dashed line gives the temperature TneLi obtained from the calculated isotope yields. 
Because of sequential feeding, it differs from the thermodynamical temperature, the un- 
corrected temperature T He Li,o being somewhat lower. Here, in order to permit the direct 
comparison with the experimental data in one figure, we display Tn e Li which has been 
corrected in the same way with the factor 1.2 suggested by the quantum statistical model. 
The model Tn e Li exhibits a more continuous rise with decreasing Zbound than the thermo- 
dynamical temperature and is in very good agreement with the measured values. We thus 
find that, with the parameters needed to reproduce the observed charge partitions, this 
temperature-sensitive observable is well reproduced. In the bin Z bound < 10 the model 
values fall below the data. Here, the experimental uncertainty is rather large but the con- 
straint on the excitation energy provided by the charge partitions is also rather weak (see 
above). The discrepancy may therefore indicate that the excitation energies for breakups 
corresponding to this bin of Z bound may be higher than assumed in the calculations. 
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The correction applied to the calculated T He Li (dashed line) results in a good overall 
agreement with the equilibrium temperature of the model calculations (full line). This 
means that the side-feeding corrections suggested by the statistical multifragmentation 
model and by the quantum statistical model are qualitatively very similar. The differ- 
ence between them is given by the deviations of the two curves. They range from about 
+20% at large Z bound to -10% at small Z hound and thus stay within the range of model 
uncertainties quoted above. They may reflect finite size effects, ignored in the quantum 
statistical model, and their variation as the spectator mass changes as a function of Z bound 
jn|. They can, equally well, be interpreted as indicating a variation of the breakup density 
with Z bound . In the region of small Z bound , i.e. at temperatures of about 8 to 10 MeV the 
correction factors of the quantum statistical model are lowered by 10% if the density is 
chosen to be pj p = 0.15 instead of 0.3. A decrease of the breakup density of that order, as 
Zbound decreases from near 40 to below 20, does not seem unreasonable. One may therefore 
expect that an analysis based on measured breakup densities (cf. ref. J38|) will bring T^Li 
into even better agreement with the equilibrium temperatures obtained with the statistical 
multifragmentation model. 

3.3 Energy spectra 

The kinetic energy spectra were studied for light charged particles up to 4 He. For the 
five species proton, deuteron, triton, 3 He, and 4 He, the spectra measured at Q[ ab = 150° 
and sorted according to Z bound are shown in fig. 2. Fit results, based on a single target 
spectator source, are given by the full lines. They were obtained by assuming that all the 
particle spectra have Maxwellian shapes 

dN/dE ~ \J E -V c ■ e ^ E - Vc ^ (7) 

where Vc is the effective Coulomb barrier and T is the temperature. The smearing of the 
spectra in the low-energy region due to the finite target thickness and a possible motion of 
the target source (at the most a few percent of the speed of light, see ref. ||) were ignored 
since their effect on the temperature parameters is too small as to be important for the 
following discussion. 

The inverse slope parameters T s i ope obtained from the fits are shown in fig. 10 as 
a function of Z bound . These kinetic temperatures also increase with decreasing Z bound , 
consistent with increasing energy deposition, but their absolute values are much higher 
than the breakup temperatures deduced from the isotope yield ratios or from the model 
description. There is also a tendency to saturate in the Z bound < 20 region. Evidently, the 
inverse slope parameters are not closely related to the breakup temperatures. 

The slope parameters of the kinetic-energy spectra calculated with the statistical mul- 
tifragmentation model are given by the lines shown in fig. 10. For the range of Z bound < 30 
they are in the vicinity of 10 MeV and nearly the same for the five particle species. This 
value is considerably higher than the internal temperatures at breakup (cf. fig. 7, full line) 
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which reflects the additional fluctuations due to Coulomb repulsion and secondary decays 
after breakup, to the extent that these effects are incorporated in the model. Apparently, 
they account for only part of the difference to the experimental slopes. Towards larger 
values of Zbound, the model results start to spread out over the range 5 to 15 MeV, with 
the protons exhibiting the lowest and 3 He the highest inverse slope parameters. 

In fig. 11, the measured kinetic energy spectra, integrated over finite ranges of Z bound , 
are shown in comparison to the model results. The two sets of experimental and model 
spectra are each normalized separately, and one overall normalization factor is used to relate 
the two sets. It was adapted to the yields of Z = 2 fragments because their calculated 
multiplicities, as a function of Zbound, are found to satisfactorily reproduce the experimental 
multiplicities reported in ref. ||. 

The main trend apparent from the comparison is a systematically increasing devia- 
tion of the experimental from the model spectra with decreasing Zbound, i.e. increasing 
centrality, and with decreasing particle mass. It not only affects the slope parameters 
describing the shape of the spectra but also the integrated intensities. The yields of hydro- 
gen isotopes, and in particular of the protons, are grossly underestimated by the statistical 
multifragmentation model. In the case of 4 He, on the other hand, the equilibrium descrip- 
tion accounts rather well for the multiplicities and kinetic energies. A major contribution 
to the observed 4 He yields is expected to come from evaporation by large fragments and 
excited residue-like nuclei which, apparently, is modelled well. 

Conceivable mechanisms that cannot explain the observed deviations include collective 
flow and Coulomb effects which both should act in proportion to the mass or charge of the 
emitted particle, contrary to what is observed. On the other hand, the commonly adopted 
scenario of freeze-out after expansion involves a pre-breakup phase during which the system 
cools not only by adiabatic expansion but also by the emission of light particles, predomi- 



nantly nucleons but also light complex particles M, H4L 45]. The spectra should reflect the 



higher temperatures at the earlier stages of the reaction, prior to the final breakup into 
fragments. In addition, there may be contributions from the primary dynamical stage of 
the reaction, not as much from the mid-rapidity source itself than from secondary scatter- 
ings with spectator nucleons. This picture is in line with the fact that higher excitation 
energies were deduced from the experimental data |9|, |13| than were needed as input for 
the calculations (fig. 8). Further investigations of the pre-breakup emission will be needed 
in order to clarify this problem (see also refs. |7], |4(| ) . 

A significant component of pre-breakup emission in the light particle yields has two con- 
sequences that deserve particular attention. The pre-breakup yields of protons, deuterons, 
and tritons are included in the double ratios used to determine the temperatures T Hep d 
and Tnedt- This violates the requirement of thermal and chemical equilibrium, which is the 
basic assumption of the method, and thus may shed doubt on the meaning of the consis- 
tency exhibited in fig. 6. On the other hand, the deduced temperatures reflect mainly the 
sensitivity of the 3 He/ 4 He yield ratio. The p/d ratio varies rather slowly with temperature, 
as evident from fig. 4 and also predicted by the quantum statistical model. Therefore, the 
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overall p/d ratio and the deduced temperatures should not be strongly affected by con- 
tributions to the hydrogen yields from earlier reaction stages. The second point concerns 
the excitation energy carried away prior to the equilibrium breakup. Part of this energy 
may be included in a calorimetric measurement of the spectator excitation. This has to be 
taken into account in the interpretation of the resulting experimental caloric curve. 

4 Conclusion and outlook 

Breakup temperatures Tn C Li, ^Hcpd; and Xfiedt were measured for target spectators in 197 Au 
+ 197 Au collisions at 1000 MeV per nucleon. In these reactions multifragmentation is the 
dominant decay channel of the produced spectator systems over a wide range of excitation 
energy and mass. The corrections for sequential feeding of the ground-state yields, based 
on calculations with the quantum statistical model, resulted in mutually consistent values 
for the three temperature observables, except in the range of very small Zbound- 

With decreasing Zb OU nd, the obtained temperatures increase from T = 4 MeV for pe- 
ripheral collisions to about 10 MeV for the most central collisions. Within the errors, the 
values for Tn e Li are in good agreement with those measured with the ALADIN spectrome- 
ter for projectile spectators in the same reaction at 600 MeV per nucleon. This invariance 
of the breakup temperature with the bombarding energy is consistent with the observed 
Zbound scaling of the mean fragment multiplicities and charge correlations and supports the 
statistical interpretation of the multi-fragment decay of highly excited spectator nuclei. 

The comparison with the results of calculations within the statistical multifragmenta- 
tion model shows that a good simultaneous agreement for the charge partitions and for the 
breakup temperatures can be achieved. A necessary requirement for a consistent statisti- 
cal description of the spectator fragmentation is thus fulfilled. The obtained equilibrium 
temperature increases less steeply with increasing centrality and stays close to T = 6 MeV 
over a wide range of Z bound , coinciding with the maximum multiplicity of intermediate- 
mass fragments. In the bin of smallest Zbound, the strong rise of the experimental Tn e Li 
to 12 MeV is not followed by the model prediction but neither is the excitation energy, 
used as input for the model calculations, well constrained for very small Zb OU nd- The side- 
feeding corrections obtained with the quantum statistical model and with the statistical 
multifragmentation model are in qualitative agreement. The remaining differences may 
largely disappear if a variation of the breakup density with impact parameter is considered 
in applying the quantum statistical model. Tn C Li for small Zbound will then be lowered 
accordingly. 

The T s i ope parameters characterizing the calculated particle spectra, although higher 
than the breakup temperatures, are still considerably smaller than the experimental values. 
The systematic increase of the deviations with decreasing particle mass indicates that 
they may be caused by light-particle emission prior to the final breakup stage. A more 
quantitative understanding of the role of the pre-breakup processes will be essential for 
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the interpretation of temperatures obtained from light-particle yields as well as of the 
excitation energies obtained from calorimetric measurements of the spectator source. 
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Figures 




Pid 



Fig. 1. Spectrum of the particle-identification variable Pid in the range 2 < Z < 4 for 
the telescope positioned at 9[ ab = 150°. The identified isotopes are indicated. 
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Energy (MeV) 

Fig. 2. Energy spectra for hydrogen and helium isotopes, measured at 9^ = 150°, for 
10-unit-wide bins in Zb oun d- All spectra are normalized with respect to each other, note 
however the scaling factors of powers of 10. The solid lines are the results of thermal- source 
fits, the dashed line indicates the trigger threshold of the telescopes for helium ions. 
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Fig. 3. Energy spectra of 6 ' 7 Li ions, measured at 9i a b = 150° (left panel), and their yield 
ratios as a function of the applied energy threshold (right panel). The energy spectra 
are integrated over the full range of Zbound, the ratios are given for four bins of Zb OU nd 
as indicated (note the linear offsets). The Maxwellian fit of the 7 Li spectrum is shown 
overlaid over the 6 Li spectrum, after an adjustment of the overall yield factor (full lines). 
The trigger threshold of the 300-/xm detector is indicated by the dashed line. 
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input 



(MeV) 



Fig. 5. Apparent temperatures TH e Li,o, 2Hepd,o, and TH e dt,o, according to the quantum 
statistical model, as a function of the input temperature Ti nput . A breakup density p/po = 
0.3 is assumed. The dotted line represents the linear relation T = T input /1.2. 
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bound 



Fig. 6. Temperatures T He Li, ^Hepd, and T Hc dt as a function of Z bound , averaged over 
bins of 10-units width. Corrections have been applied as described in the text. The 
error bars represent the statistical uncertainty. The systematic uncertainty, caused by the 
extrapolation of the yields of helium isotopes below the identification threshold, is indicated 
by the brackets. For clarity, T Hc Li and T He dt are laterally displaced by 2 units of Z bound . 
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Fig. 7. Temperatures Tn e Li of the target spectator from the present experiment at E/A 
= 1000 MeV (dots) and of the projectile spectator at E/A = 600 MeV (open squares) as a 
function of Abound- The data symbols represent averages over bins of 10-units width and, for 
clarity, are laterally displaced by 1 unit of Zb oun d- Statistical and systematic contributions 
are included in the displayed errors. The lines are smoothed fit curves describing the 
breakup temperature T breakup (full line) and the isotopic temperature T He Li (dashed line) 
calculated with the statistical multifragmentation model. Note that the trigger threshold 
affected the data of ref. |J at Zb OU nd > 65. 
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Fig. 8. Excitation energy E x /A as a function of the mass A for the ensemble of excited 
spectator nuclei used as input for the calculations with the statistical multifragmentation 
model. The area of the squares is proportional to the intensity. 
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Fig. 9. Mean multiplicity of intermediate-mass fragments {Mimf) (top) and mean charge 
asymmetry (ai 2 ) (bottom) as a function of Zb OU nd, as obtained from the calculations with 
the statistical multifragmentation model (open circles) in comparison to the experiment 
(dots, from ref. ||). The dashed and dotted lines show the results of the calculations with 
excitation energies E x /A 15% above and 15% below the adopted values, respectively. Note 
that the trigger threshold affected the data of ref. at Zbound > 65. 
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bound 

Fig. 10. Inverse slope parameters T s i ope for hydrogen and helium isotopes, obtained from 
fits to the measured (data symbols, 9i ab = 150°) and calculated (lines) energy spectra, as 
a function of Z bound . 
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Fig. 11. Energy spectra, measured at 9i a b = 150°, of light charged particles p, d, t, 3 He, 
and 4 He for three intervals of Zb oun d as indicated. The dots represent the measured spectra, 
the open circles are the results of the calculations with the statistical multifragmentation 
model. The spectra are normalized as stated in the text. 
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